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Compact stars with strong magnetic fields (magnetars) have been observationally determined to 
have surface magnetic fields of order of 10 14 — 10 15 G, the implied internal field strength being several 
orders larger. We study the equation of state and composition of dense hypernuclear matter in strong 
magnetic fields in a range expected in the interiors of magnetars. Within the non-linear Boguta- 
Bodmer-Walecka model we find that the magnetic field has sizable influence on the properties of 
matter for central magnetic field B > 10 17 G, in particular the matter properties become anisotropic. 
Moreover, for the central fields B CI > 10 19 G, the magnetized hypernuclear matter becomes unstable, 
which limits the range of admissible fields in magnetars to fields below the critical value B CT . 

I. INTRODUCTION 

Anomalous X-ray pulsars and soft y-ray repeaters are observationally identified with highly magnetized neutron 
stars with surface magnetic field ~ 10 14 — 10 15 G [H-Q . This class of compact stars has been conjectured theoretically as 
"magnetars" 043; *- e -; neutron stars which posses magnetic fields that are by several orders of magnitude larger than 
the canonical surface dipole magnetic fields B ~ 10 12 -10 13 G of the bulk of pulsar population, which are commonly 
deduced from the magnctic-dipolc radiation model of pulsar spin-down in combination with measured spin and spin- 
down rates. The integral properties of magnetars, i.e., mass, radius, moment of inertia, etc will depend sensitively 
on the equation of state of matter in strong magnetic fields, if the central fields are sufficiently strong. Furthermore, 
other processes, such as the cooling and the magnetic field evolution will sensitively depend on the composition of 
matter in strong magnetic fields. Fermionic matter in strong magnetic field experiences two well-known quantum 
mechanical phenomena: the Pauli paramagnetism and the Landau diamagnetism. The first is due to the interaction 
of the spin of the fermion with the magnetic field and, therefore, is relevant for both charged and uncharged fermions. 
The second phenomenon is relevant only for charged particles, and is particularly strong for light particles, which in 
the case of compact stars are the leptons. 

Neutron star sequences may feature massive objects with masses M ~ 2M Q . These massive compact stars are 
likely to develop cores which are composed of matter that differs from the ordinary nuclear matter composed of only 
neutrons and protons. One possibility is that heavy baryons (hyperons) will appear once the Fermi energy of neutrons 
becomes of the order of their rest mass. Although the hyperons were considered even before the discovery of pulsars 
and their identification with the neutron stars |8(, their presence in the cores of neutron stars is still quite uncertain, 
for different theoretical models predict quite different outcomes for the maximum mass of hypernuclear compact stars, 
often in contradiction with the known empirical data on pulsar masses. The models that are based on rclativistic 
density functional methods [9l-[l2| predict masses that are not much larger than the canonical mass of a neutron star 
which clearly contradicts modern observations. Masses of the order of < I.SMq were obtained in non-relativistic 
phenomenological models (l3l Ibij , while microscopic models based on hyperon-nucleon potentials, which include the 
repulsive three-body forces, predict low maximal masses for hypernuclear stars [l5l.[l6j. The avenues for reconciliation 
of the large pulsar mass and the hyperonization (and more generally strangeness) of dense matter have been explored 
recently in Refs (HUH. 

The influence of strong magnetic fields on the highly dense electron gas in the context of neutron star matter was 
studied in Rcfs. |29l - [33i . The strong magnetic field effects on dense nuclear matter (n,p, e system) have been studied 
previously in Refs. |34| - |40| . In Refs. [4l|, [42j the effect of strong magnetic field on hypernuclear matter has been 
studied. The structure of strongly magnetized neutron star branch of compact objects was studied for different field 
configurations (toroidal, poloidal, etc) in Refs. (43l - |45j . In the case of the adjacent branch of white dwarfs, strong 
magnetic fields were found to lead to highly super-Chandrasekhar mass (M ~ 2.3 — 2.6Mq) white dwarfs (46l - |48j . 
which can be related to the observed features of a number of peculiar Type la supernovae [49( ■ 

The purpose of this paper is a unified study of the properties of dense nuclear matter which feature both hypernuclear 
degrees of freedom and strong magnetic fields. We will focus here on the equation of state of matter and its composition. 
In doing so, we will concentrate on the effect of Landau quantization of the motion of fermions in strong magnetic 
fields, which modifies the single particle energies of charged fermions. Furthermore, as well known, the magnetic 
field contributes directly to the energy density and the pressure of matter, which is limited by the magneto-static 
equilibrium of the stars. In the simplest form it can be formulated for uniform self-gravitating fluids [50(. In the 
case of neutron stars the Chandrasekhar- Fermi limiting field strength is ~ 10 18 G [5l[. Similar limits were obtained 
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for stable dense nuclear and quark matter, based on ferromagnetic symmetry breaking [52j . However, the uniform 
field distribution, while giving a crude order of magnitude estimates, certainly requires refinements, since the density 
profiles of neutron stars are not uniform. Consequently, the local magneto-static equilibrium will be achieved for 
different field strength and any dynamo mechanism generating fields will carry the imprint of inhomogeneous density 
profile in the star. Furthermore, the processes of supernova collapse will leave behind a strongly non-uniform field 
distribution of frozen-in fields. Thus, more realistic treatment of the equation of state of matter requires inclusion of 
some physically motivated field profiles of the stars such as magnetars. (We note that local equation of state is not 
affected by the density profile of the field, and therefore depends only on the local value of the magnetic field strength 
and induction, see Ref. (13). Below we shall consider a more realistic situation in which the profile of the magnetic 
field is density dependent, where the latter is modelled by some reasonable functional form with proper limiting value 
at the surface of the star consistent with the observations. 

This paper is organized as follows. In Sec. [TT] we discuss our model for hypernuclear matter in strong magnetic 
fields. Our numerical results are presented in Sec. IIIII Section ITVl contains a summary of our results. 

II. MODEL 

As discussed in the introduction, the treatment of hypernuclear matter within nuclear models which are based 
on different principles has lead to array of results, some of which are inconsistent with measurements of pulsar 
masses. Relativistic models of hypernuclear matter are flexible enough to be adjusted to the current phenomenology, 
therefore they provide a suitable framework for extensions which incorporate the influence of strong magnetic fields. In 
particular, we consider the non-linear Boguta-Bodmer-Walecka model which is based on relativistic density functional 
theory [jjl |53T - |57| . This model is able to describe ground state phenomenology of nuclear matter and nuclei (see, e.g., 
\E§L ) • We anticipate that the main conclusions and results of the present study will not change qualitatively if a 
different model of dense hypernuclear matter is chosen. 

The Lagrangian density of hypernuclear matter in a static magnetic field (hereafter also i?-ficld) is given by 

C s = C m + £ f , (1) 
where the matter part of the Lagrangian is given by 

£m = ^ i>B {ilftD^ - m B + g aB o- - g^Bl^ - 9 p b1^t b • p M ) ipB 
B 

+\d^a - \mW - U(a) 

--^^ v + -mlw^ - l -p^ v ■ pT + \m 2 p p p ■ p^ 

+ ^2^i(i-r^-mi)^ l , (2) 

where the index B labels baryons present in the matter, ipB,' l Pi,o~,uj and p arc fields of baryons, leptons, cr-mesons, 
w-mesons and p-mesons, with masses m.B , ?nj , , m w and m p respectively; g a B,gu>B and g p s are coupling constants 
for interactions of a, ui and p mesons respectively with the baryon B. The scalar self interaction term in the matter 
Lagrangian is 0, [UH^ 

U(cr) = i gi m N {g aN a) s + ^ 92 (ffo-jver) 4 , (3) 



and 



u MU = dyuj^ - d^ v , (4) 
/V = dvPp - dnp v , (5) 
D» = d» +ieQA». (6) 



We choose the field vector potential gauge as = (0, — yB, 0, 0), where B is the magnitude of magnetic field and eQ 
the charge of the particle, e the (positive) unit of charge. For this particular gauge choice the magnetic field is along 
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the z axis of Cartesian coordinate system, i.e., B = Bz. In the mean field approximation, the baryons acquire the 
effective masses 

m* B = m B - g aB cr, (7) 
where a is given by its ground state expectation value 



< 



1 / v-^ ( B) dU \ . 

= ^ a \2^9.Bn s —J, (8) 

with the scalar density 

(-B) 

9 rPF ml 

=d 3 PB , (9) 



(2^) 3 Jo 



p% + m\ 



where ps is the momentum and p F the Fermi momentum of the baryon B. The electromagnetic field Lagrangian 
density is given by 

C t = ~F iaf F^ 1 (10) 

where is the electromagnetic field tensor. 

The total energy density and pressure of the system can be obtained by computing the energy-momentum tensor 
from the Lagrangian density (fTJ). The result is 

T" v = T£" + Tf , (11) 



-f 

where the matter part is given by 

T% = e rn u^u v - P(g^ - u»u v ) + hM^ x F^ + AI ,yX F^), (12) 

with e rn being the matter energy density, P - the thermodynamic pressure, - the magnetization tensor. The 

field part of the energy-momentum tensor is given by 

T ^ = _l_ F ^ F U + ±_gl i V F PC F (13) 

J 47T 1D7T 

In the following we will neglect the electric field, as there are no macroscopic charges in the bulk matter. Therefore, 
Eqs. (fT2|) and (fT3|) reduce, respectively, to [60L rolj 



m 




rpflV 

1 f 


4tt V 



B»B V 



(14) 



- -A - (15) 
2 J J 4tt ' K ' 

with B^B^ = —B 2 and M being the magnetization per unit volume. 

In the presence of the magnetic field, the motion of the charged particles is Landau quantized in the direction 
perpendicular to the magnetic field. We choose the coordinate axes in such a way that B is along z-axis. Then, the 
single particle energy of any charged particle at n-th Landau level is given by 



E n = ^pl + m 2 + 2ne\Q\B, (16) 

where m is the mass of the particle, p z is the component of momentum along z direction. The Landau levels in 
Eq. (fPoj) assume integer values n = 0,1,2... for spin-up states and n = 1,2,3... for spin-down states for positively 
charged particles. For negatively charged particles n takes on values n = 0,1,2,... for the spin-down states and 
n = 1, 2, 3... for spin-up states. The zero-temperature number density of charged baryons and leptons is given by 

n O = E ( 2 " S nfi )^/ P 2 F -2ne\Q\B, (17) 
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with 



2e\Q\B J ' 



(18) 



where pp is the Fermi momentum. At zero temperature the number density of neutral baryons is expressed via their 
Fermi momentum as 



n N 



P F 

3tt 2 



(19) 



The chemical potentials of baryons are given by fis = P*b + ^'fJuB + P3I3 9pB, where 



Mb = yp^'+m**, (20) 
with m* B being the effective mass of the baryon B in the mean field approximation. We will further need the quantities 

1 11 z * 



3 n B, 



(21) 
(22) 



P B 



( B) 

which are the ground state expectation values of u> and ps fields with 1^ being the isospin projection for the baryon 
B. The chemical potentials of leptons are given by 



Mi 



P^p 2 + TO 2 



(23) 



where p F ^ and mi are the Fermi momenta and masses of leptons, respectively. Now we are in the position to write 
down the matter energy density defined by the matter Lagrangian density, C m , given by Eq. @ 



Stt 

e\Q\B 

(2tt) 2 

e\Q\B 



B ^mi»* B -mi ln P ^ )+/ ^ 



rt mai 

x; E(2-<w>) 



in 



B 



B' n=0 



PB>{n)n* B , + (m* B , + 2ne\Q\B) In 



pi(n)fii + (mf + 2ne\Q\B) In 



PB>{n) + ii* B , 



m* B , + 2ne\Q\B , 



l—e,/i n—0 



pi{n) + pi 



^mf + 2ne\Q\B 



(24) 



where the sums over i? and £?' include the uncharged and charged baryons, respectively, and the following short-hand 
notation has been used 



V{n) = sJp 2 F -2ne\Q\B. 
The matter pressure can be easily constructed from Eq. ([24| using the thermodynamic relation 

P = E MB??B + E Mi"/ - £rn- 



(25) 



(26) 



We now account for the fact that the matter in compact stars is charge neutral and in beta equilibrium. The first 
requirement relates the partial densities of charged particles according to 



n p + n E + = n s - + n s - + n e + n,, 



(27) 
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Equilibrium with respect to weak interactions further implies the following relations 



Mp = M« - Me, Mm = Me, MA = Mn, ME" = Mn + Me, 

Ms = Mn, Ms+ = Mp, Me- = M» + Me, Ms = Mn- 
These conditions are normalized such that the total number of baryons is reproduced, i.e. 



(28) 



(29) 



The solutions of the field equations at any given baryon density ns and zero temperature are found in the mean-field 
approximation under the constraints (|27[) . (|28p . and (|29[) . Subsequently, the energy density and the thermodynamic 
pressure of the matter are obtained using Eqs. (f2~4")) and (|2l)|) . 

In the rest frame of the hadronic fluid, with B field along the z axis, the matter and field parts of the energy-stress 
tensor are given, respectively, by 



( e m 
P -MB 




V 



P-MB 



(30) 







P 



1 f 




(31) 



The total energy density of the system is given by the sum of the matter and field contributions 

It is seen from Eqs. ([50)) and (|3"Tj) that the the pressure in the perpendicular direction to the magnetic field is 

B 2 

P ± =P-MB+—, 
Sir 

and the pressure in the direction parallel to the magnetic field is given by 

B 2 

P\\=P • 

11 8tt 



(32) 



(33) 



(34) 



The couplings in the hypernuclear Lagrangian are fixed as follows. In the nucleonic sector the nucleon-meson coupling 
constants are chosen according to Refs. 0, I551457} and reproduce the bulk properties of nuclear matter (the binding 
energy E/B = —16.3 McV, the saturation density uq — 0.153 fm -3 , the asymmetry energy coefficient a asy = 32.5 
MeV, the incompressibility K = 240 MeV and the effective nucleon mass at the saturation m* jm = 0.8). In the 
hyperonic sector, the hyperon-w coupling constants is derived from the SU(6) symmetry of the quark model (6214641 ] . 
To infer the hyperon-a coupling constants we write the potential depth of a hyperon Y in nuclear matter at saturation 
as 



and assign certain value to the potential Uy, which then fixes the value 



(35) 

We take the potential depth for A-hyperon 
Ua = —30 McV, as obtained from the analysis of A-hypernuclei [62j, [65j. For S-hyperon we assign Us = —18 MeV 
according to experimental data of Refs. [HH, [66| ■ According to the recent E-hypernuclei data the potential for E 
hyperons is repulsive (6?J and we adopt the value Us = 30 McV. 

To model the density profile of the magnetic field we adopt the following formula [68| 



B 



B x +BA1 



exp 



(36) 



As discussed in the introduction the present formula allows us to model the realistic physical situations where the 
magnetic field of the star is non- uniform which is a realistic situation. The functional dependence of the profile (|36[) on 




FIG. 1: Left panel: Variation of total pressure as a function of normalized baryon number density for fixed magnetic fields 
B c = and B c = 10 18 G and several field profiles, /3 = fO" 3 , 7 = 6 (dashed lines), (3 = f(T\ 7 = 4 (dashed-double-dotted 
lines), and /3 — > 00, i.e., B c = constant (dashed-dotted lines). For each pair of curves the upper branch is for P± and the lower 
branch for P11. Right panel: Same as in the left panel, but for fixed j3 = 10 _1 , 7 = 1 and for two values of B-field, B c = 10 18 G 
(dashed lines) and B c — 3 x 10 18 G (dashed-dotted lines). The solid line corresponds to the case B c — 0. 



density is constructed such as to take into account the fact that the magnetic field in the cores of compact stars may 
be higher than that at the surface. The parameters (3 and 7 control the relaxation from the central value B c to the 
asymptotic value at the suface B s . Obviously, the parameter /3 controls the amount of field decay at the saturation 
density, whereas the parameter 7 controls the width of the transition. Obscrvationally inferred value of the surface 
magnetic field of magnetars is of the order of 10 15 G. It is plausible to asume that the central field could be in the 
range of 10 17 — 10 18 G. Therefore, we will consider a number of different values of B c , while keeping B s = 10 15 G. 
We have verified that the results are insensitive to the precise value of the surface field B s . 



III. RESULTS AND DISCUSSION 



Our model of the equation of state of hypernuclear matter in strong magnetic field is parameterized in terms of four 
physical quantities: the boundary values of the field at the center and the surface, the amount of its decline from the 
center to the saturation density (which is close to the crust-core interface) and the steepness of the transition between 
the two asymptotic values of the B-field. The latter two are described by the parameters j3 and 7, respectively. While 
the surface £>-ficld has been fixed at the value of 10 15 G, we have explored a range of central magnetic fields; it turns 
out that the equation of state is unaffected by the magnetic fields below the value B c < 10 18 G. 

Figure [T] shows the equation of state of hypernuclear matter in strong magnetic field. The left panel shows the 
dependence of pressure on the density for a fixed central 5-field and various field profiles. For non-zero magnetic 
field the pressure splits into the parallel and transverse components, which is evident from Eqs. and and is 
due to the gauge field contribution. The case with fi = 10 _1 acquires a larger field at the saturation density n = n , 
therefore the transition to the asymptotic value of B c occurs earlier than in the case of f3 = 10 -3 , both cases having 
approximately the same width of transition (see Fig. [2]). It is clearly seen that the low-density behavior of the equation 
of state with constant £>-ficld implies unrcalistically large anisotropy of magnetic field up to the surface of the star, 
which is inconsistent with the inferred surface magnetic field of magnetars B s ~ 10 15 G. The right panel of Fig. Q] 
shows the equation of state in case of fixed density profile /3 = 10 _1 and 7 = 1, but for various central magnetic fields. 
The main effect seen in the figure is the splitting of the parallel and transverse pressure as the field becomes larger, 
this effect being sizable for fields above 10 18 G. Furthermore, we see that the increase of B c causes the profile for P11 
to become softer; the fact that P11 < P±_ is due to negative contribution of the field pressure to P|| . 

From the left panel of Fig. [T] we observe that there is an onset of matter instability for chosen combinations of /? 
and 7 (corresponding to different magnetic field profiles) for fixed central field B c . This instability is induced in the 
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FIG. 2: Upper panels: Dependence of Pu and B on the normalized baryon number density for different magnetic field profiles 
and B c = I0 18 G. The dots show the reference case B c — 0. The solid and dashed lines correspond to = 0.1 and O.OOf, 
respectively. For each /3 we choose a pair of 7's; in the first case we have 7 = 1 and 7 = 4, whereas in the second case 7 = f 
and 7 = 6. The first 7 value for each f3 corresponds to a less steeply rising curve. Lower panels: Same as in the upper pannel 
but for B c — 3 x 10 18 G and 7 = 1 and 7 = 1.4 for /3 = 0.1 and 7=1 and 7 = 3.5 for /3 = I0 -3 . The specific values of 7's were 
chosen to demonstrate the onset of instability. 



P11 component of the pressure, which decreases with the density, because its field profile implies an increase of field 
strength. 

To study the onset of instability in a more systematic manner, consider the dependence of P|| on density for two 
central values of magnetic field B c = 10 18 G and B c = 3 x 10 18 G and several sets of 7 and (3 parameters. This is 
shown in the upper and lower panels of Fig. [2] It is seen that for a given value of /3 the corresponding EoS becomes 
softer as the 7 is increased. Consequently, beyond a certain critical value of 7 and in a certain density regime Pu ceases 
to increase and subsequently decreases with the further increase in rn,. This implies that matter becomes unstable 
above that value of density for that particular B c and magnetic field profile. This is evident from Fig. [3J where also 
the magnetic field profiles corresponding to the onset of instability are shown. For comparison we also show results 
for each f3 with the minimum value of 7 taken to be 1. Note that the maximum value of 7 is taken such that Pu forms 
a plateau as a function nj,. Furthermore, it is evident that with the decrease of (3, the instability occurs at larger 
values of 7 and rib, which is also evident from Fig [2] 

The instability arises due to the negative contribution from the field energy density (pressure) to the pressures of 
magnetized baryons and leptons in the direction to the magnetic field, which is evident from Eq. (|34p . Since for any 
particular B c and magnetic profile, the field strength increases with the increase of rib, more negative contribution is 
added to Py with the increase of n&. Consequently, at a certain density, Py ceases to increase and then decreases with 
the increase of nt. The critical values of 7 for the onset of instability as a function of log /3 are shown in Fig. [3J The 
potential instability discussed above was found in dense baryonic matter in Ref. [H] and has been confirmed since 
by Ferrer et al. [70, [zj and Dexheimer et al. [72[ . Similar behavior was also found in strange quark matter [6l], [z3 ■ 
Occurrence of instability in nuclear matter in the presence of large magnetic field constrains the value of central 
magnetic field of the magnetars. 

To demonstrate the effects of Landau quantization on the EoS we show in Fig. [4] the EoS as a function of density 
with the field energy density B 2 /8ir subtracted. It is seen that the effects of Landau quantization are insignificant 
even at field values 3 x 10 18 G and become sizable only for fields of order 10 19 G. 

Here we should mention that the true minimum which determines the ground state of the matter is controlled by 
the Gibbs potential, which contains an additional term M B. However, even at the strong field strength considered 
in the present study, the magnetization effect is negligible. For example, the contribution of the magnetization to Pj_ 
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FIG. 3: Dependence of y c - the c 
the field B c = 3 x 10 18 G, the da; 




FIG. 4: Dependence of matter pressure on normalized baryon number density without the B 2 /8ir term, which demonstrates 
the effect of Landau quantization. 



in Eq. (|33[) is at least an order of magnitude smaller than that due to matter and field, depending upon the strength 
of the field and density. Hence, in our numerical results, we do not consider the effect of magnetization. However, 
it should be noted from Eqs. (|33|) and (|34| that, in principle, the magnetization leads to anisotropy of the total 
pressure (apart from that due to the magnetic field itself) of the system and contributes negatively to the pressure in 
the perpendicular direction of the magnetic field. Therefore, inclusion of magnetization further strengthens the main 
conclusion of this study. 

In this context, it is also worth mentioning that the effect of interaction of the anomalous magnetic moment (AMM) 
of particles with field could stiffen EoS [42|. It could, in principle, dominate over the softening effects by Landau 
quantization for fields strength above 5 x 10 18 G i.e., well above those considered in this paper. It will start dominating 
over the Landau quantization effect for fields ~ 5 x 10 19 G. However, at such a high field strength the field energy 
density contribution is more than enough to overwhelm the AMM interaction. Furthermore, the effect of AMM 
interaction is density dependent and is most significant at low densities. In this respect, it is important to take 
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FIG. 5: Abundances of different species as functions of normalized baryon number density without (left panel) and with (right 
panel) magnetic field. 



a realistic profile of a star, which will suppress the contribution from AMM interactions because it implies low fields 
at low densities. 

Finally, it is interesting to ask how the strong i?-fields affect the composition of matter. In Fig. [5] we show the 
composition of hypcrnuclear matter as functions of matter density in cases with and without B-field. It is seen that 
the B-field has an effect only on the lightest particles, i.e., the leptons. For these we see "Landau oscillations" in 
their population as a function of the field. These oscillations reflect the population of new Landau levels with the 
increase of the density. 

IV. SUMMARY 

We have studied the influence of strong magnetic fields on the high density nuclear matter, including the possibile 
appearance of hyperons at large densities. In doing so we have implemented realistic density profiles of magnetic 
fields, which assume that the field decreases from the center of the star and reaches asymptotically its value at the 
surface ~ 10 15 G characteristic for magnetars. The functional dependence of the field is taken as in Ref. [f38j]. We have 
found that the magnetic field effect is significant when B c > 10 17 G and at large field strength the matter becomes 
unstable. This is due to the negative contribution of the field pressure to the baryonic and leptonic pressures in the 
direction parallel to the magnetic field, which renders the total pressure of the system anisotropic. We have found 
that the onset of instability depends on the magnetic field profile (parameterized in terms of /3 and 7) as well as on 
the central field value B c . The instability depends on particular values and form of the density profile, but it sets in 
always for critical central field values B CI « 10 19 G for any values of (3 and 7. This gives a natural bound for the 
central magnetic field of neutron stars. 

Acknowledgement 

M. S. acknowledges the support by the Alexander von Humboldt Foundation. B. M. acknowledges partial support 
through research Grant No. ISRO/RES/2/367/10-11. We thank Dr. X. G. Huang for discussions. 



[1] G. Vasisht and E. V. Gotthelf, Astrophys. J. 486, L129 (1997), arXiv:astro-ph/9706058. 

[2] C. Kouveliotou, S. Dieters, T. Strohmayer, J. van Paradijs, G. J. Fishman, C. A. Meegan, K. Hurley, J. Kommers, I. Smith, 

D. Frail, et al., Nature (London) 393, 235 (1998). 
[3] P. M. Woods, C. Kouveliotou, J. van Paradijs, K. Hurley, R. M. Kippen, M. H. Finger, M. S. Briggs, S. Dieters, and G. J. 

Fishman, Astrophys. J. 519, L139 (1999), arXiv:astro-ph/9903267. 



10 



[4] R. C. Duncan and C. Thompson, Astrophys. J. 392, L9 (1992). 
[5] V. V. Usov, Nature (London) 357, 472 (1992). 

[6] C. Thompson and R. C. Duncan, Mon. Not. Roy. Astron. Soc. 275, 255 (1995). 

[7] C. Thompson and R. C. Duncan, Astrophys. J. 473, 322 (1996). 

[8] V. A. Ambartsumyan and G. S. Saakyan, Astron. Zh. 37, 193 (1960). 

[9] N. K. Glendenning, Astrophys. J. 293, 470 (1985). 
[10] N. K. Glendenning and S. A. Moszkowski, Physical Review Letters 67, 2414 (1991). 
[11] F. Weber, ed., Pulsars as astrophysical laboratories for nuclear and particle physics (1999). 
[12] A. Sedrakian, Progress in Particle and Nuclear Physics 58, 168 (2007), arXiv:nucl-th/0601086. 
[13] S. Balberg, I. Lichtenstadt, and G. B. Cook, Astrophys. J. 121, 515 (1999), arXiv:astro-ph/9810361. 
[14] H. Djapo, B.-J. Schaefer, and J. Wambach, Phys. Rev. C 81, 035803 (2010), 0811.2939. 
[15] M. Baldo, G. F. Burgio, and H. . Schulze, ArXiv Astrophysics e-prints (2003), arXiv:astro-ph/0312446. 
[16] I. Vidafia, D. Logoteta, C. Providencia, A. Polls, and I. Bombaci, EPL (Europhysics Letters) 94, 11002 (2011), 1006.5660. 
[17] L. Bonanno and A. Sedrakian, A&A 539, A16 (2012), 1108.0559. 
[18] W.-Z. Jiang, B.-A. Li, and L.-W. Chen, Astrophys. J. 756, 56 (2012), 1207.1686. 
[19] I. Bednarek, P. Haensel, J. L. Zdunik, M. Bejger, and R. Mahka, A&A 543, A157 (2012), 1111.6942. 
[20] M. Oertel, A. F. Fantina, and J. Novak, Phys. Rev. C 85, 055806 (2012), 1202.2679. 
[21] K. Masuda, T. Hatsuda, and T. Takatsuka, ArXiv e-prints (2012), 1205.3621. 
[22] R. Lastowiecki, D. Blaschke, H. Grigorian, and S. Typel, ArXiv e-prints (2011), 1112.6430. 
[23] S. Schramm, R. Negreiros, J. Steinheimer, T. Schiirhoff, and V. Dexheimer, ArXiv e-prints (2011), 1112.1853. 
[24] V. Dexheimer, J. Steinheimer, R. Negreiros, and S. Schramm, ArXiv e-prints (2012), 1206.3086. 
[25] S. Weissenborn, D. Chatterjee, and J. Schaffner-Bielich, Nuclear Physics A 881, 62 (2012), 1111.6049. 
[26] S. Weissenborn, D. Chatterjee, and J. Schaffner-Bielich, Phys. Rev. C 85, 065802 (2012), 1112.0234. 
[27] F. Gulminelli, A. R. Raduta, J. Margueron, P. Papakonstantinou, and M. Oertel, ArXiv e-prints (2012), 1209.0270. 
[28] R. Mallick, ArXiv e-prints (2012), 1207.4872. 
[29] V. Canuto and J. Ventura, Fundam. Cosmic Phys. 2, 203 (1977). 
[30] I. Fushiki, E. H. Gudmundsson, and C. J. Pethick, Astrophys. J. 342, 958 (1989). 
[31] A. M. Abrahams and S. L. Shapiro, Astrophys. J. 374, 652 (1991). 

[32] I. Fushiki, E. H. Gudmundsson, C. J. Pethick, and J. Yngvason, Annals of Physics 216, 29 (1992). 

[33] O. E. Roegnvaldsson, I. Fushiki, E. H. Gudmundsson, C. J. Pethick, and J. Yngvason, Astrophys. J. 416, 276 (1993). 

[34] S. Chakrabarty, D. Bandyopadhyay, and S. Pal, Physical Review Letters 78, 2898 (1997), arXiv:astro-ph/9703034. 

[35] D. Bandyopadhyay, S. Chakrabarty, P. Dey, and S. Pal, Phys. Rev. D 58, 121301 (1998), arXiv:astro-ph/9804145. 

[36] A. Broderick, M. Prakash, and J. M. Lattimer, Astrophys. J. 537, 351 (2000), arXiv:astro-ph/0001537. 

[37] I.-S. Suh and G. J. Mathews, Astrophys. J. 546, 1126 (2001), arXiv:astro-ph/9912301. 

[38] A. K. Harding and D. Lai, Reports on Progress in Physics 69, 2631 (2006), arXiv:astro-ph/0606674. 

[39] W. Chen, P.-Q. Zhang, and L.-G. Liu, Modern Physics Letters A 22, 623 (2007), arXiv:astro-ph/0505113. 

[40] A. Rabhi, C. Providencia, and J. Da Providencia, Journal of Physics G Nuclear Physics 35, 125201 (2008), 0810.3390. 

[41] P. Yue, F. Yang, and H. Shen, Phys. Rev. C 79, 025803 (2009), 0902.3699. 

[42] A. E. Broderick, M. Prakash, and J. M. Lattimer, Physics Letters B 531, 167 (2002), arXiv:astro-ph/0111516. 

[43] M. Bocquet, S. Bonazzola, E. Gourgoulhon, and J. Novak, Astron. Astrophysics 301, 757 (1995), arXiv:gr-qc/9503044. 

[44] C. Y. Cardall, M. Prakash, and J. M. Lattimer, Astrophys. J. 554, 322 (2001), arXiv:astro-ph/0011148. 

[45] J. Frieben and L. Rezzolla, ArXiv e-prints (2012), 1207.4035. 

[46] U. Das and B. Mukhopadhyay, Phys. Rev. D 86, 042001 (2012), 1204.1262. 

[47] U. Das and B. Mukhopadhyay, ArXiv e-prints (2012), 1205.3160. 

[48] A. Kundu and B. Mukhopadhyay, Modern Physics Letters A 27, 1250084 (2012), 1204.1463. 

[49] R. A. Scalzo, G. Aldering, P. Antilogus, S. Aragon, C. Bailey, C. Baltay, S. Bongard, M. Buton, C. Childress, N. Chotard, 

Y. Copin, H. K. Fakhouri, et al., Astrophys. J. 713, 1073 (2010), 1003.2217. 
[50] S. Chandrasekhar and E. Fermi, Astrophys. J. 118, 116 (1953). 
[51] D. Lai and S. L. Shapiro, Astrophys. J. 383, 745 (1991). 

[52] A. Perez Martinez, H. Perez Rojas, and H. J. Mosquera Cuesta, International Journal of Modern Physics D 17, 2107 

(2008), 0711.0975. 
[53] J. D. Walecka, Annals of Physics 83, 491 (1974). 
[54] J. Boguta and A. R. Bodmer, Nuclear Physics A 292, 413 (1977). 
[55] N. K. Glendenning, Physics Letters B 114, 392 (1982). 

[56] N. K. Glendenning, Zeitschrift fur Physik A Hadrons and Nuclei 326, 57 (1987). 

[57] N. K. Glendenning, Zeitschrift fur Physik A Hadrons and Nuclei 327, 295 (1987). 

[58] S. A. Chin, Annals of Physics 108, 301 (1977). 

[59] B. D. Serot and J. D. Walecka, Adv.Nucl.Phys. 16, 1 (1986). 

[60] V. R. Khalilov, Phys. Rev. D 65, 056001 (2002). 

[61] X.-G. Huang, M. Huang, D. H. Rischke, and A. Sedrakian, Phys. Rev. D 81, 045015 (2010), 0910.3633. 
[62] C. Dover and A. Gal, Comments Nucl.Part.Phys. 12, 155 (1984). 

[63] J. Schaffner, C. B. Dover, A. Gal, C. Greiner, D. J. Millener, and H. Stocker, Annals of Physics 235, 35 (1994). 
[64] J. Schaffner and I. N. Mishustin, Phys. Rev. C 53, 1416 (1996), arXiv:nucl-th/9506011. 

[65] T. Fukuda, A. Higashi, Y. Matsuyama, C. Nagoshi, J. Nakano, M. Sekimoto, P. Tlusty, J. K. Aim, H. En'yo, H. Funahashi, 



11 



et al., Phys. Rev. C 58, 1306 (1998). 
[66] P. Khaustov, D. E. Alburger, P. D. Barnes, B. Bassalleck, A. R. Berdoz, A. Biglan, T. Burger, D. S. Carman, R. E. Chrien, 

C. A. Davis, et al., Phys. Rev. C 61, 054603 (2000), arXiv:nucl-ex/9912007. 
[67] S. Bart, R. E. Chrien, W. A. Franklin, T. Fukuda, R. S. Hayano, K. Hicks, E. V. Hungerford, R. Michael, T. Miyachi, 

T. Nagae, et al., Physical Review Letters 83, 5238 (1999). 
[68] D. Bandyopadhyay, S. Chakrabarty, and S. Pal, Physical Review Letters 79, 2176 (1997), arXiv:astro-ph/9703066. 
[69] M. Sinha and B. Mukhopadhyay, ArXiv e-prints (2010), 1005.4995. 

[70] E. J. Ferrer, V. de La Incera, J. P. Keith, I. Portillo, and P. L. Springsteen, Phys. Rev. C 82, 065802 (2010), 1009.3521. 

[71] L. Paulucci, E. J. Ferrer, V. de La Incera, and J. E. Horvath, Phys. Rev. D 83, 043009 (2011), 1010.3041. 

[72] V. Dexheimer, R. Negreiros, and S. Schramm, ArXiv e-prints (2011), 1108.4479. 

[73] X.-G. Huang, A. Sedrakian, and D. H. Rischke, Annals of Physics 326, 3075 (2011), 1108.0602. 



